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We construct the integration measure over the moduli space of an arbitrary number of N kinds 
of dyons of the pure SU(N) gauge theory at finite temperatures. The ensemble of dyons governed 
by the measure is mathematically described by a (supersymmetric) quantum field theory that is 
exactly solvable and is remarkable for a number of striking features: 1) The free energy has the 
minimum corresponding to the zero average Polyakov line, as expected in the confining phase; 
2) The correlation function of two Polyakov lines exhibits a linear potential between static quarks 
in any iV-ality non-zero representation, with a calculable string tension roughly independent of 
temperature; 3) The average spatial Wilson loop falls off exponentially with its area and the same 
string tension; 4) At a critical temperature the ensemble of dyons rearranges and de-confines; 5) The 
estimated ratio of the critical temperature to the square root of the string tension is in excellent 
agreement with the lattice data. 
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I. INTRODUCTION 



Isolated dyons in the pure Yang-Mills theory are (anti) self-dual solutions of the equation of motion D^F^ = 0, 
which in an appropriate gauge are Abelian at large distances from the centers and carry unity electric and magnetic 
charges with respect to the Cartan generators C m : 

1 r 

zh E m = B m = — - — — C m - (1) 

I |r| 

For the SU(N) gauge group on which we focus in this paper there are N — 1 Cartan generators associated with the 
simple roots of the group, C m = diag(0, 1, —1, 0, 0) where the 1 is on the 771 th place, that are supplemented by 
the N th generator Cn — diag(— 1, 0, 1) to make the set of N dyons electric- and magnetic-neutral. The first N — 1 
dyons are also called the Bogomolny-Prasad-Sommerfield (BPS) monopoles The last N th dyon is sometimes 
called the Kaluza-Klein monopole: in the gauge where the first N—l dyons are described by a static field the last 
has time-dependent fields inside the core. However, by a periodic time-dependent gauge transformation one can make 
the last one time-independent (at the cost of a time dependence inside other dyons) , therefore the distinction of the 
last dyon is illusory and we shall treat all N of them on the same footing. The cyclic symmetry of N dyons is evident 
from the D-branes point of view 0, Q . 

In this paper we explore the properties of a semiclassical vacuum built of a large number of dyons of N kinds. 

To make the semiclassical calculation of the Yang-Mills partition function well defined one needs 1) to expand 
about a true saddle point of the action, ii) to be sure that the quantum fluctuation determinant is infrared-finite. 
For an isolated dyon the first is true but the second is false. For an arbitrary superposition of N different-kind dyons 
the second is true but the first is wrong. To satisfy both requirements one can consider N dyons as constituents of 
the Kraan-van Baal-Lee-Lu (KvBLL) instantons with non-trivial holonomy 0, @, which are saddle points of the 
Yang-Mills partition function as they are exact self-dual solutions of the equations of motion. At the same time the 
fluctuation determinant about the KvBLL instanton is infrared-finite (and actually exactly calculable Q) since its 
total electric and magnetic charges are zero. 

KvBLL instantons generalize the standard Belavin-Polyakov-Schwartz-Tuypkin (BPST) instantons [7J having triv- 
ial holonomy. The mere notion of dyons and of the KvBLL instantons (also called calorons) alike imply that the 
Yang-Mills field is periodic in the Euclidean time direction, as it is in the case of non-zero temperature. Therefore 
we shall be considering the Yang-Mills partition function at finite temperature. However, the circumference of the 
compactified space can be gradually put to infinity, corresponding to the zero temperature limit. In that limit, the 
temperature can be considered as an infrared regulator of the theory needed to distinguish between the trivial and 
non-trivial holonomy. After all the temperature of the Universe is 2.7 K ^ 0. 

In this context, the holonomy is defined as the set {/i m } of the gauge-invariant eigenvalues of the Polyakov loop L 
winding in the compactified time direction, at spatial infinity: 



L = P exp 1 1 j 1 dt A t \ = V diag (e 2 ^ 1 , e 2wi ^ , . . . , e 2 ^ N ) V' 1 , ^ /x m = 0. (2) 

By making a global gauge rotation one can always order the eigenvalues such that 

Mi < M2 < • • • < MiV < MiV+i = Mi + !■ (3) 
which we shall assume done. If all eigenvalues are equal up to an integer, implying 

triv ( 4r — 1 when TO < k, t r 1 Ar I a\ 

mv = 1 n - > where A; = 1,...,N, (4) 

I W w]i en 777 > k 

the Polyakov line belongs to the SU(N) group center, and the holonomy is then said to be "trivial": L tl i V — 
diag (exp(2Trik / N) , ...,exp(2ivik/N)) , k = 1, .... N. Standard BPST instantons, as well as their genuine periodic gen- 
eralization to non-zero temperatures by Harrington and Shepard [8( , possess trivial holonomy, whereas for the KvBLL 
instantons the gauge-invariant eigenvalues of the Polyakov line assume, generally, non-equal values corresponding to 
a "non-trivial" holonomy. Among these, there is a special set of equidistant /i m 's that can be named a "maximally 
non-trivial" holonomy, 

^^ri^- (5) 

having a distinguished property that it leads to Tr L — 0. Since the average Polyakov line is zero in the confining 
phase, the set ([5]) can be also called the "confining" holonomy. 
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Whatever the set of /it's is equal to, it is a global characterization of the Yang-Mills system. Integrating over all 
possible p,'s is equivalent to requesting that the total color charge of the system is zero @. Since ju's are constants, 
the ultimate partition function has to be extensive in these quantities, meaning Z = exp[— F({/j,}, T)V] where F is 
the free energy and V is the 3-volume. If f ({//}) has a minimum for a particular set {fi}, integration over /x's is 
done by saddle point method justified in the thermodynamic limit V — ► oo, hence the Yang-Mills system settles at 
the minimum of the free energy as function of /x's. The big question is whether the pure Yang-Mills theory has the 
minimum of the free energy at the "confining" holonomy ([5]), or elsewhere. 

It has been argued long ago [9( that the pure Yang-Mills theory has N minima of the free energy at the trivial 
holonomy (|4]). To that end, one refers to the perturbative potential energy as function of spatially constant A4 (9l.[l0|: 

(r, \2rp3 JL 

= V (Mm ~ Vn?[l ~ (Urn " Mn)] 2 , (6) 

which, indeed, has N minima (all with zero energy) at the trivial holonomy ([3]). corresponding to N elements of the 
center of SU(N). The "confining" holonomy ([5]) corresponds, on the contrary, to the non-degenerate maximum of 
([6]), equal to 

ppert, m ax_ T , (27r) 2 T 3 iV 4 -l 

180 TV 2 1 ' 

The large volume factor in eq. ([6]) seemingly prohibits any configurations with non-trivial holonomy, dyons and KvBLL 
instantons included. 

A loophole in this dyon- killing argument has been noticed in Ref. If one takes an ensemble of dyons, with their 
number proportional to the volume, the non-perturbative dyon-induced potential energy is also proportional to the 
volume and may hence override the perturbative one, possibly leading to another minimum of the full free energy. 
This scenario was made probable in Ref. [6( where it was shown that the non-perturbative potential energy induced 
by a dilute gas of the KvBLL instantons prevailed over the perturbative one at temperatures below some critical T c 
estimated through A, the Yang-Mills scale parameter, and that the trivial holonomy was not the minimum of the full 
free energy anymore. Below that critical temperature the KvBLL instantons dissociate into individual dyons. The 
problem, therefore, is to build the partition function for dissociated dyons, and to check if the full free energy has a 
minimum at the confining holonomy ([5]). We get an affirmative answer to this question. 

The moduli space of a single KvBLL instanton is characterized by 4iV parameters (of which the classical action 
is independent); these can be conveniently chosen as 3c? coordinates of N dyons constituting the instanton, and 
their U(l) phases, 3N + N = AN. In the part of the moduli space where all dyons are well separated, the KvBLL 
instanton becomes a sum of N types of BPS monopoles with a time-independent action density. At small separations 
between dyons the action density of the KvBLL instanton is time-dependent and resembles that of the standard 
BPST instanton. The KvBLL instanton reduces to the standard BPST instanton in the two limiting cases: i) trivial 
holonomy (all /u's are equal up to an integer) and any temperature, ii) non-trivial holonomy but vanishing temperature, 
provided the separations between dyons shrink to zero as ~ p 2 T where p is the standard instanton size. 

The quantum weight of the KvBLL instanton is determined by a product of two factors: i) the determinant of 
the moduli space metric, ii) the small-oscillation determinant over non-zero modes in the KvBLL background. The 
latter has been computed exactly in Ref. [f| for the SU(2) gauge group; recently the result has been generalized to 
any SU(N) The former is also known exactly (see the references and discussion in the next Section). These 

achievements, however, are limited to the case of a single KvBLL instanton with unity topological charge. To build 
the dyon vacuum, one needs multi-instanton solutions, with the topological charge proportional to the volume, similar 
to the case of the instanton liquid model [ill, EH ■ Although there has been serious progress recently in constructing 
general multi-instanton solutions with non-trivial holonomy and their moduli space metric [3], a desirable explicit 
construction is still lacking. 

Nevertheless, the moduli space metric can be constructed by combining the metric known for N different-kind 
dyons of the SU(N) group with another known metric for same-kind dyons, and by taking into due account the 
permutational symmetry between identical dyons. One of the two ingredients, the metric for different-kind dyons, 
is known exactly for all separations and involves only Coulomb-like interactions. The other ingredient related to 
the same-kind dyons is more complex. The metric for any separations between same-kind dyons allows for charge 
exchange between dyons and involves elliptic functions of separation. For two dyons of same charge the exact metric 
was found by Atiyah and Hitchin (l5j from the requirement that the Riemann tensor constructed from the metric 
must be self-dual. When the separation between same-kind dyons exceeds their core sizes, charge exchange dies out 
exponentially with the separation, the metric becomes simple and can be written for any number of same-kind dyons 
with the exponential precision [la , Il7| . It involves only the Coulomb-like interactions and is in fact very similar to 
that for different-kind dyons but with opposite signs in the Coulomb bonds. It is this marriage of the asymptotic 



4 



form of the metric for same-kind dyons, valid with exponential precision, with the metric for different-kind dyons, 
valid for any separations, that we shall explore in this paper. 

In fact, it may prove to be sufficient for an accurate description of the ensemble of dyons in the thermodynamic 
limit (V — > oo), as dyons of the same kind repulse each other whereas dyons of different kind attract each other. 
Therefore small separations between same-kind dyons, where our measure is only approximate, may be statistically 
unimportant. [Unfortunately, taking the large N limit does not help since only nearest neighbors in color interact, 
hence at any N > 2 there are only twice more different-kind bonds than same kind.] Indeed, we find that, despite an 
approximate integration measure, the Lorentz symmetry is, in a sense, restored at T — > 0: the "electric" string tension 
as determined from the correlation of Polyakov lines turns out to be independent of temperature and equal to the 
"magnetic" string tension determined from the area law for spatial Wilson loops, for all representations considered. 
The free energy itself also has a reasonable limit at T — > 0. However, to remain on the safe side, we claim the results 
only for sufficiently high temperatures (but below the deconfinement phase transition) where dyons do not overlap on 
the average, and the metric used is justified. 

This study is exploratory as we ignore many essential ingredients of the full Yang-Mills theory. In particular, we 
consider only the ensemble of dyons of one duality, whereas CP invariance of the vacuum requires that there must 
be an equal number of self-dual and anti-self-dual configurations, up to the thermodynamic fluctuations ~ \f~V [l3| . 
We basically ignore the determinant over non-zero modes, taking from it only certain known salient features, like 
the renormalization of the coupling constant and the perturbative potential energy Our aim is to demonstrate 
that the integration measure over dyons has a drastic, probably a decisive effect on the ensemble of dyons, that 
the ensemble can be mathematically described by an exactly solvable field theory in three dimensions, and that the 
resulting vacuum built of dyons has certain features expected for the confining pure Yang-Mills theory. 

A more phenomenological and lattice-oriented attempt to describe the ensemble of the KvBLL instantons, has been 
proposed recently in Ref. [l~8] . 

II. INTEGRATION MEASURE OVER DYONS 
A. Different-kind dyons 

The metric of the moduli space of a single SU(N) KvBLL instanton, written in terms of N different-kind dyons' 
coordinates and U(l) pha ses, was first conjectured by Lee. Weinberg and Yi generalizing the previous work by 
Gibbons and Manton [17] . and then confirmed by Kraan [20] by an explicit calculation of the zero- mode Jacobian 
from the Atiyah-Drinfeld-Hitchin-Manin-Nahm (ADHMN) construction [HI, HI] for the SU(N) caloron [H]. It was 
later checked independently in Ref. [24J , also by an explicit calculation of the Jacobian. 

There are several equivalent ways to present the metric of a single SU (N) KvBLL instanton; we use here the form 
suggested originally by Gibbons and Manton [Tt} (although these authors considered another case - that of same kind 
dyons, see the next subsection): 

ds 2 = G + (dip m + W mm > ■ dxm>) G m ^ {dipn + W nn > ■ dxn>), m,n — 1...N. (8) 

Here x m , m = 1...N, are the 3d centers of dyons, and ip m are their U(l) phases. G mn is a symmetric matrix composed 
of the Coulomb interactions between dyons that are nearest "neighbors in color" : 

f~i I / A . i 1 1 I dm.n — 1 ^m,n+l / n x 
Gmn = O mn I 47W m + + I - - -, (9) 

\ P^m X m _ 1| X-m+ll/ X m _j_i| ^m— 1| 

where v m = /i m +i — ji m \ V\ + ... + = 1 (see Section I). Periodicity in color indices is implied throughout the paper: 
m = N + 1 is equivalent to m = 1, m = is equivalent to m — N. W m „ are three N x N symmetric matrices 
composed of the electric charge-magnetic charge interaction potential w(x): 

W m „ = 5 mn (w(x m - x ro _i) + w(x m - x m+ i)) - 5 TO) „_iw(x m - x m+ i) - 5 m ,„+iw(x m - x m _i) (10) 

where w(x) satisfies the equation e abc <9;,w c = — x a /|x| 3 ; to solve it, one has to introduce a Dirac string singularity. 
Choosing the string, e.g., along the third axis and parameterizing the separation vector between two dyons in spherical 
coordinates x = |x|(sin0 cos<^, sin# sin^), cos#) one finds w(x) = (— cot# sin</>, cot# coscj), 0)/|x|, such that the 
factor (dip + w • dx) written for relative coordinates combines into 

dS 3 = dip + cos 6 d(f> (11) 

which is a familiar 1-form encountered, e.g., in the theory of rigid-body rotations, being the third projection of the 
angular velocity. The fact that this quantity arises in the context of electric charge-magnetic charge interaction, was 
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known for quite a while [2a] . A more general, basis-independent way to present the interaction is via the Wess-Zumino 
term 



The temperature factors T have been dropped in eqs.([8] fT0|) but can be anytime restored from dimensions. We stress 
that the metric ([5]) is exact and valid for any 3d separations between dyons inside the KvBLL instanton, including 
the case when they strongly overlap. For more details on this metric see Refs. [20l. l2~4l| . 

The integration measure for the KvBLL instanton is 



N 



Y\_ d 3 yL m dip m \J det g, \J det g = det G, (12) 



m— 1 



where g is the full AN x AN metric tensor given by eq. ([8]). In computing det g one can shift d%\) m — » dip' m 
dipm + W m fc • dxfc, therefore, detg = (detG) 3 detG -1 = (detG) 2 , hence y/detg = det G where G is given by 



eq. © [24(. In that reference it was also shown that in the limit of trivial holonomy or small temperatures, the 



integration measure (fl"2")) reduces to the well-known one for the standard BPST instanton [2j 

B. Same-kind dyons 

For multi-KvBLL instantons, a new element appears, namely two or more same-kind dyons are present. For two 
dyons of the same kind, the metric splits into a flat 4c? space for center-of-mass coordinates and a non-flat Ad space M-i 
for relative coordinates r and tp. Self-duality implies that M 2 is a self-dual Einstein manifold. Gibbons and Pope [28[ 
proposed the following form for the metric: 

ds 2 = fdr 2 + a 2 d^l + b 2 d^j + c 2 d^l, y/g = \fabc\sin6, (13) 

where 

rfSi = — sin tp d9 + cos tp sin 9 d<j>, 
dTii = cos ip dO + sin tp sin 6 d<j>, 

and dT,3 is the third 'angular velocity' (fTT|) ; a,b,c,f are functions of the dyon separation r. Self-duality requires that the 
Riemann tensor built from the metric f) 13[) satisfies Ra/s^s = \ \J~§ e-ySnA R a p kX wri ich leads to the system of first-order 
equations 

\ da b 2 + c 2 — a 2 

= A, and cyclic permutations of a, b, c, (14) 

/ dr 2bc 

where A = or 1. The value A = 1 is chosen from symmetry considerations |2!|. Eqs. (|14p have a simple solution fiH ] 



2a „ / 2a 2aB , . 

/ = -/3 A /l+_, a = b = (3rJl+—, c= (15) 

with any a,/?. We fix a, (3 from the asymptotics of the metric of two dyons of the same kind m at r — > oo: 
j3 = y/2nv m , a = —l/(2irv m ) [3(1]. To get the full 8x8 metric tensor in terms of the two dyons' coordinates 

xi = X + ±r, dtpi = rf* + \dip, 
x 2 = X - |r, dijj 2 =d$>- \dtp, 

we add to eq. (|T3)) the flat metric for center-of-mass coordinates 87ri/(dX-G!X + d^ d^) and obtain the metric which 
is very similar to eq. ([8]): 

ds 2 = Gij dxi-dxj + (d^i -dxii) G^ 1 (d-ipj +Wjj> -dxj>), i,j = 1,2, (16) 
_ ^ / Airv m | Xl _ X2 | |xi-x 2 | | /i 7 \ 

G y - 2 4 2 , (J-7) 

\ |x!-x 2 | ^ Vm |xi-x 2 | / 

Wn = ( ^ w ( x ^- x ^ w(xi-x 2 ) \ w(x) = 1 ( _ . . . 

13 \ w(xi-x 2 ) -w(xi-x 2 ) r K ' x v ; v ; 
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Note the opposite sign of Coulomb interactions in eq. (fl7|) as compared to eq. ([9]). 

It is easy to generalize eqs. (ll6HrS|) to any number of dyons, all of the m th kind. One extends the summation in 
eq. (|16jl from i,j<2 to i,j < K, where K is the number of same-kind dyons, and modifies eqs. (|17H8l) as 

a f^-E^^, i = j, w |-2E^w(x 4 -x fc) , i = j, 

{ R^TT' I 2w(x i -x J ), 1 + 3 

Eqs. (|16ll9[) were derived by Gibbons and Manton [ItJ (with other coefficients related to another scale convention) 
from considering the classical equations of motion for K identical monopoles at large separations. 

Although eq. (fl~5|) from where the metric (fT6|) stems, is an exact solution of the Einstein self-duality equation, it is 
believed that eq. (|16p is valid only for large separations r > 1 /2irv m T (we restore the explicit temperature factor here) . 
Note that a very similar metric (JSJ) for different-kind dyons is proven to be valid at any separations. A somewhat 
superficial reason for the difference between same- and different-kind dyons was noted in Ref. [l9j |: while the metric 
(j8]| is positive-definite, the metric (fT6|) goes to zero at r = l/2-Kv m T. A deeper reason is that while the metric (|8|) 
describes a system with total electric and magnetic charges zero (the KvBLL instanton), the metric (fTB)) is applied 
to a system where neither is zero. 

A non-trivial solution of the Einstein self-duality eq. (I14|) was found by Atiyah and Hitchin (AH) [15J (for more 
details see [29j). In the AH solution, the a,b,c functions are given by elliptic integrals, and a(r) ^ b(r). It follows 
then from eq. (|13[) that the relative U(l) orientation angle ip enters the metric in an essential way, in particular, a 
shift of i\) is not an isometry anymore. The functions a, b, c of the AH solution differ from those of the solution (fT5|) 
by terms of the order of ~exp(— Aitv m Tr) which die out exponentially at large separations [T(ij |. Therefore, at large 
r the AH solution takes the form of eq. (JTSJ) such that ip enters the metric in a trivial way , as in eq. (fTfj)) . 

For the AH solution, the metric determinant goes to zero at even larger r — l/(4v m T) [TU, Ojj]. Physically, this point 
corresponds to an axially-symmetric two-monopole solution where two monopoles coincide. When dyons overlap, what 
should be called "separation" becomes ambiguous; r is defined only in the context of a concrete parameterization of 
the field. In the ensemble, the zero of the metric determinant means a vanishing contribution to the partition function, 
actually imposing a very strong repulsion between same- kind dyons. The same is true for the approximate metric 
(|T5| which we shall use below. One can think that the approximate metric p^|) . because of the strong repulsion it 
imposes, strongly suppresses in the statistical mechanics' sense configurations with close same-kind dyons, where the 
metric becomes inaccurate. In other words, the approximate metric (|16[) may be accurate for statistically important 
configurations. This hypothesis needs a detailed study, of course. Its consequences, however, turn out to be reasonable. 



C. Combining the metric for same-kind and different-kind dyons 



The explicit form of the metric tensor for K KvBLL instantons made of N kinds of dyons is not known (for the 
latest development see Ref. fT3|). Below we suggest an Ansatz for this metric, satisfying the known requirements. 
One has to combine the metric (|5]) for N different-kind dyons with that for K same-kind ones. The solution of the 
problem is almost obvious if one takes the approximate metric (|16l) for same-kind dyons, as it has exactly the same 
form as the metric ((HJ for different-kind dyons. Since the metric cannot 'know' to which instanton a particular dyon 
belongs to, it must be symmetric under permutations of any pair of dyons of the same kind. Important, the metric 
of the moduli space of self-dual solutions must be hyper-Kahler [2i| . 

Let indices m,n = 1...N refer to the dyon kind (or 'color') and indices k = 1...K number the dyons of the same 
kind. The coordinates of the i th dyon of the m th kind are (x mi , Vw) = Vmi = V% a = 1, 2, 3, 4. To shorten notations, 
we introduce instead of the multi-index {mi) a single index A = (mi) running from 1 to KN. 

We write the full metric tensor as 



ds z = gAa.sp dy^dy 1 ^ = G AB dx A -dx B + (dip A + W AA ,-dx A ,) G AB (dip B + W BjB / -dx B 



(20) 



where, explicitly, 



GAB — Gmi,nj — $mn$ij I 47TZ/ m + ^ 



3 m , n — 1 



+ 1 
" ^rn — 1 . 



E 



X-mk | 



(21) 
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W AB = W mit nj = S mn 5ij y^w(x mj -x m _ lifc ) +y^w(x rm --x m+1 ,fc) -2y]w(x mi -x mfe ) (22) 

^m,n— 1 w(x m j X m ^x_j) <5m,n+iw(x m j X m _ij) -)- 2J mn w(x m ^ X 7T1 j ) | ^-^-^ • 

The inverse matrix G~ AB in (I20p is understood according to the relation G~ a x c Gcb = $ab = fimndij- 

Note that the Coulomb bonds in eq. (|21|) for the same-kind dyons have an opposite sign from those for neighbor 

kind, and have a twice larger coefficient. The coefficients, —1, 2, —1, are actually the scalar products of simple roots of 

the SU(N) group, supplemented by an additional non-simple root to make the matrix cyclic-symmetric. This remark 

allows the generalization of eqs. (|2H22[) to any Lie group. 

The constructed AKN x AKN metric tensor gAa,B/3 is hyper-Kahler. It means that there exist three "complex 

structures" 1(a), a = 1, 2, 3, (all three are AKN x AKN matrices) such that 

I{a)g = gl(af ("T" means transposed) (23) 

and which satisfy the Pauli algebra, 

I(a)I(b) = e abc I(c) - S ab l. (24) 
Related to 1(a), there are three Kahler symplectic 2- forms 

w(o) = r>(a)mc 7 dy% A dy^, il(a) = -0(a) T , (25) 

where 

0(a) = I{a)g. (26) 



The 2-forms ui(a) are closed: 



dw(a) = or - j ^Sl{a) B p < c 1 dylAdy B Adyl = Q. (27) 
°Va 



Explicitly, the three Kahler forms to (a) have the same form as in Ref. [17| for same- kind dyons, only Gab and W AB 
should be now taken from eas. (|2H22|) : 

Lu(a) = 2(dip A + W AA ,-dx A ,) Adx a A - G BC e abc dx b B A dx c c . (28) 



With G AB and given by eas. (|21I22|) . the three Kahler forms u>(a) (or Vl(a) B p y c~i) are fixed from eq. 

and the complex structures 1(a) are found from inverting eq. (f2"6")l . We have checked that the algebra (|23l24p is then 
satisfied for any choice of w(x) in eq. (|22p. It is the closure of the two- forms, eq. ([27)) , that requests that w(x) is the 
electric charge-magnetic charge interaction potential satisfying the equation e abc db~w c = — x a /|x| 3 . 

We note further properties of the constructed G AB and ~W AB : 

• symmetry: G AB = G B a, W A b = W Bj4 , meaning, of course, G mijnj = G n ^ mi , W mi , nj = W n]tmt 

• overall "neutrality": J2nj G mi,nj = 47w m , J2mi G mi,nj = ^v n , V (( , W ,,,,.„, = 0, J2mi W mi,nj = 

• identity loss: dyons of the same kind are indistinguishable, meaning mathematically that det G is symmetric 
under permutation of any pair of dyons (i<-+j) of the same kind m 

• factorization: in the geometry when dyons fall into K well separated clusters of TV dyons of all kinds in each, 
det G factorizes into a product of exact integration measures for K KvBLL instantons, det G = (det G\) K where 
G\ is given by eq. ([9]). 

The integration measure over the moduli space of K KvBLL instantons of the SU (N) gauge group is 

K N 

IT IT / dXmi d ^ mi V det 9, y/detg = det G. (29) 

i— 1 m — 1 

In deriving the last relation we notice that in the determinant, d?p m i can be shifted by *W m i,m'i' ■dx. m 'i', hence 
det g = (detG) 3 detG -1 = (detG) 2 , therefore Vdet g = detG where the KN x KN matrix G is given by eq. ((21) . 
We have also checked this result by an explicit calculation of the determinant of the full AKN x AKN metric tensor 
g. Since G is independent of the U(l) angles ip m i, integration over tp can be omitted. 
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D. Dyons' fugacity 



Fugacity is a term from statistical mechanics of grand canonical ensembles (where the number of particles is not 
fixed) denoting the weight with which a particle contributes to the grand partition function. Let there be K m dyons 
of the m th kind, m = 1...N. For a neutral system of K KvBLL instantons the number of dyons of every kind is equal, 
K\ = ... = Kn = K, however we shall consider the general case of non-equal K's for the time being: one can always 
project to the neutrality condition. For an arbitrary set of AT's, G is a (Kx + ... + K^) x (K\ + ... + K^) matrix given 
by eq. (2}. 

We write the partition function of the grand canonical ensemble as a sum over all numbers of dyons of each kind: 



N K. 



Z = E K j l K i iill/ ( d3 *™f) detG(x), (30) 

where / is the x-independent factor - the fugacity - accompanying every integral over x. Since det G(x) is symmetric 
under permutation of same- kind dyons, the identity factorials are needed to avoid counting the same configuration 
more than once. If one likes to impose the overall neutrality condition, viz. that only configurations with the equal 
number of dyons of different kind contribute to the partition function (K\ = ... = A/v = K) one integrates eq. (|30[) 
over auxiliary angles: 

1 X II II /(rf 3 x™/)detG(x) (31) 



^ / n 2tt "' 2tt K x \ K n \ 

K 1 ...K N J ° 1 JV m=l i=l 



N K 



Et^Iv II Ufa*«f) detG(x). 

K ^ ' m=l i=l J 



We shall see below, however, that the neutrality condition will be taken care of dynamically, therefore the additional 
integration pip is in fact unnecessary. 

As fugacity is x-independent it can be established from the limit when all dyons are grouped into A-plets of 
different-kind dyons, forming infinitely dilute neutral KvBLL instantons, such that the measure factorizes into a 
product of individual instanton measures. The latter is known to be [2~j ] 



N 
m— 1 



d 3 x m detGi(x)2 2JV 7r 3A ' (32) 



where Gi(x) is the N x N matrix ([9]) for one KvBLL instanton, see eq. IT2"]) . This must be multiplied by the factor 

(/x 4 /T) j \J2-Kg 2 coming from 4./V zero modes of the instanton. Here n is the ultraviolet cutoff and g is the 

bare coupling constant given at that cutoff @, |24j. Multiplication by this factor makes eq. (|32p dimensionless, as it 
should be. In addition, eq. (|32p is multiplied by the exponent of minus the classical action of the instanton, equal to 
(A/fi)~ N where A is the Yang-Mills scale parameterizing the coupling constant in the Pauli-Villars regularization 
scheme, and by the dimensionless factor [det(— -D 2 )] where D 2 is the Laplace operator in the instanton background. 
The last factor arises from integration over non-zero modes; it is understood that the small-oscillation determinant 
is normalized to the free (zero field) determinant and UV-regularized by the Pauli-Villars method. It is known that 

the normalized and regularized [det(— -D 2 )] 1 is proportional to , times the exponent of minus the perturbative 
potential energy Q , times a slowly varying function of dyon separations [1, [3l| . 

Combining all factors we observe that the Pauli-Villars mass \x cancels out (as it should be in a renormalizable 
theory) and we obtain the dyon fugacity 



f=Y~g~* C (33) 



where c is proportional to A~3 ; it is made dimensionless by a combination of temperature and dyon separations. The 
relative (un)importancc of c in the dynamics of the ensemble, as compared to the measure factor (f3"0")) , is illustrated by 
the powers of A: their ratio is (—1/3) : 4 = —1/12. For the time being we shall set c — 1 and recall it in the discussion 
in Section IV. The coupling g 2 in (f3"3"l) starts to 'run' at the two-loop level not included here. Ultimately, its precise 
argument is determined self-consistently from the action density of the ensemble [32|- In the study of the large- A'' 
behavior it will be important that c = 0(1) whereas 1/g 4 = 0(N 2 ), hence the fugacity / = 0(N 2 ). 
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III. DYON PARTITION FUNCTION AS A QUANTUM FIELD THEORY 

We now face an interesting problem of finding the correlation functions in the ensemble of dyons whose grand 
partition function is given by eq. (|30[) . The renormalized Yang-Mills scale parameter A creeps in via the fugacity 
([55]) . therefore all physical quantities will be henceforth expressed through A. The temperature also enters explicitly 
via eq. (|33|) ; the temperature factors are understood in all Coulomb bonds in the matrix G (j21j) as well, to make them 
dimensionless. Thus, the partition function and the ensuing correlation functions depend, generally, both on A and 
T. 

The ensemble defined by the partition function (|30[) is a very unusual one, as it is governed by the determinant of a 
matrix G whose dimension is equal to the number of particles, and not by the exponent of the interaction energy, as is 
common in statistical mechanics. Of course, one can always write detG = expTrlogG = exp(L r ; n t), however then the 
interaction potential U\ n t will contain not only 2-body, but also 3,4,5... -body forces that are increasingly important. 
At the same time, the statistical mechanics of an ensemble governed by the determinant-induced interactions can be 
transformed into an equivalent quantum field theory which considerably simplifies its handling. 

To that end, we first notice that a matrix determinant can be presented as a result of the integration over a finite 
number of anticommuting Grassmann variables [33| , 

det G = [ J] difa d^ A exp U\ G AB 1> B ) , (34) 

where the usual convention [33| for anticommuting integration variables is understood: 

ipAipB + 4>BipA = 0, ip A i/> B + Wb^A = °> Va^b + 4>b*I> a = for any A, B, 

dtl>\ d^ A = 0, / d^ A dip A $ A i> A = 1. (35) 



In our case A = (mi) is a multi-index where m — X...N is the dyon kind, and i = l...K m is the number of a dyon of 
the m th kind. We rewrite identically the partition function (|3"0)) as 



N K 



Z= H K ; K i II II / d3x ™ / #m, exp G m ^ nJ l/} nj J (36) 

where G m i, n j is a matrix made of Coulomb interactions, eq. (f2Tj) , and / is the fugacity (|33[) where we put c = 1. 
Having obtained G in the exponent, it is now possible to express its Coulomb matrix elements from path integrals, 
extending the Polyakov trick [34j to anticommuting variables. 

A. Off-diagonal elements: ghost fields 

We first present the off-diagonal (i ^ j) elements of ip^Gip by means of a functional integration over anticommuting 
(or ghost) fields. In the next subsection we present the diagonal (i — j) elements with the help of a functional 
integration over commuting (boson) fields. 

Let us consider the Gaussian path integral over N anticommuting fields Xm (x) coupled to the anticommuting source 
J2i bPmi5(x--Xmi) ~ -iAm-M^(x-Xm-M)]: 

N 

y = I] / DxlDxm exp fdxJ2 { 7- dxl-dXm (37) 

+ [(^miXm(x) + xL( x )^ ml ) <5(X-X ml ) - (vL-l,iXm(x) + xLW^m-l,i) ^(x-X m _i,i) I . 



Although we do not write it explicitly to save space, we assume that y is normalized to the same path integral with 
the kinetic term but without the source term. The subscript m is periodic: m — N + 1 is equivalent to m = 1 and 
m = means m = N. 
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The path integration of an action that is quadratic in anticommuting variables is performed in the same way as the 
Gaussian path integral over bosonic variables, with the result 



y = exp 



mirm-l : j 



v mi A mj | 



X 



771—1, J | 



mj | 



| ^777—1 , i , j | 



(38) 



Owing to the cyclic summation over m, the last term doubles the first one, and we correctly reproduce the off-diagonal 
(* 7^ j) P ar t of i^ mi G m i^ n ji\) n j in eq. p6|l (cf. the second line in eq. (|2Tj) ). However, the sum in eq. (|38p contains an 
extra diagonal divergent term 2tp mi ip m i/\x. m i — x. m i\ which is absent in eq. (|36[) and hence should be canceled. 

B. Diagonal elements: boson fields 

Next, we present the diagonal (i—j) part of tp^Gip by means of a Gaussian integration over bosonic fields v m , w m . 
Let us consider 



• dv m ■ dw r , 



X = IT / Dv mDw m exp d^y^ l 

+ ^ ^V'LV'rm<5(x-X mi ) - ^^m-litx-X^ulj B ro (x) + (J(x-X mi ) - ^X-^i)) W m (x) 



(39) 



To make this path integral formally convergent one assumes that the integration over either v m or w m goes along 
the imaginary axis. As in the case of ghost fields in the previous subsection, we do not write it explicitly but assume 
that X is normalized to the same path integral with the kinetic term but without the source term (the second line in 
eq. (|39]>). 

Integrating (|3"9"|) over w m we obtain a functional 5-function: 



whose solution is 



^m(x) 



|X X m _X 5 i 



^ "m(x) = 0. 



(40) 



The Jacobian following from the 5-function, det (— (T/47r)<9 2 ) , cancels with the same Jacobian from the normalization 
integral. Substituting v m (x) from eq. (|40p back into eq. (|39[) and using the cyclic symmetry of the summation over 
m we obtain 



X = exp 



| Xmi X m _i j | | X m _|_i j | 



(41) 



and 



The divergent term at i — j, namely — 2f/r ni ?/> m j/|x m j— x m j|, cancels exactly the unwanted extra term in eq. 

we reproduce precisely the diagonal term in ip mi G m i tn ji r 'nj of eq. (|36|) (cf. the first line in eq. (|2T|) K 

Thus, we have fully reproduced the factor exp(?/^ Gip) in the partition function (|36p with the help of the integration 
over anticommuting ghost (xtuXm) and ordinary boson (v m ,w m ) variables. The Coulomb interactions have been 
traded for kinetic energy terms of those fields. Apparently, v m are [/(l)-^ -1 Abelian electric potentials, and w m are 
their dual fields. 
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C. Synthesis: the equivalent quantum field theory 
We now use eqs. (|37l39[) to rewrite identically the partition function (j3"6]) . We have 

T 



Z = Y[ j D xln D Xm Dv m Dw m exp J 



m—l 

00 pKi 



-in 



(dxL'dXm + dv m -dw m ) 



x y - 

Ki=0 



f_ 



dxi / dipldipi exp 4twiV4^i + #1 (xi(^i)-X2(^i)) + i (x\ (xi)~xU x i)) V'l 



+ («i(a; 1 )-U2(a;i)) ViV'i + (wi(zi)-w 2 (a:i)) 



(42) 



A„=0 



fA A 



dxiy j dip N dip n exp 



Anvffipjsr'ipN + {xn(x n )-xi(x n )) + i (xn(x n )-x\(xn)) 



+ {vn(xn)-vx(x n )) ipljipN + (w n (x n )-w-l(x n )) 



In writing eq. (|42p we have used that K m identical integrals over dx mi dip mi dip m i appear in the partition function, 
where i = l...K m is a 'dumb' index labeling integration variables. Therefore one representative of such integral for 
every dyon kind m is taken to the power K m . 

In each line in eq. (|42|) integration over dip^dipm can be trivially performed, given the rules (|35|) : it reduces to 
expanding the exponents in eq. (|42|) to the terms bilinear in ip^, ip m . For example, for m — l we get 



J^\\ f J dx [ AlTVl + Ui-X2)(xi-X2)(x) + («i - «a)(x) 

47ri/i + (Xi-Xa)(Xi-X2) + Oi - V2) 



Ai=0 

= ex P f / J dx 
and similarly for other values of m. We obtain 

N 



Ki 



Z = f[ j D xln D Xrn Dv m Dw m CXp J I — {dXm ' d Xm + Qv m ■ dw r , 

m—l m ^ 

}■ 



(43) 



(44) 



+ / UnV m + ( v rn - V m+ i) + (xL - xL+l)(Xm ~ Xm+l) 



Given the cyclic symmetry in the summation over m, the last line can be rewritten in a more nice way. We introduce 
the function 



N 



,W m —Wn 



(45) 



and recall that v m = /i m +i — [L m where [i m are the eigenvalues of the Polyakov line, see the Introduction. The second 
line in eq. (|4"2"|) can be written as 



/ 



dT d 2 T 

(-47T/Z m + V m )- h Xln T. q Xn 

dw m dw m dw n 



(summation over repeated indices is understood) where 



dw m 
d 2 T 

dWmdWn 



W m -W m + 1 _ W m -i-W„ 



W rn -W m + 1 _|_ W m -1 — W, 
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The final result for the dyon partition function is 



Z = jDx^ DxDvDw exp Jd 3 xl^^ (dixlndiXrn+diVrndtWm) 



(46) 



/ 



dT d 2 T 



Eq. (|46p should be divided by the normalization integral being the same expression but with zero fugacity /. In fact 
the normalization integral is unity and can be omitted. Indeed, integrating over v m gives 5(— (T/4ir)d 2 w m ) whose 
only solution is w m = const, whereas the Jacobian is det _1 (— {T/An)d 2 ). This Jacobian, however, is immediately 
canceled by the integral over the ghost fields Xm- Therefore, the quantum field theory defined by eq. (|46|) is the full 
result for the dyon partition function. 



IV. GROUND STATE: 'CONFINING' HOLONOMY PREFERRED 

The fields v m enter the partition function (|46[) only linearly. Therefore, they can be integrated out right away, 
giving rise to a J-function 

Dv m — ► S[-^d 2 w m + f^-). (47) 



This S- function restricts possible fields w m over which one still has to integrate in eq. (|46p . Let w m be a solution to 
the argument of the (5-function. Integrating over small fluctuations about w gives the Jacobian 



T „, „ . d 2 T 



Jac = dct 1 -— d 2 8 mn + f 



47r dw m dw n 



(48) 



Remarkably, exactly the same functional determinant (but in the numerator) arises from integrating over the ghost 
fields, in the same background w: 



Jdx'Dx exp Jd 3 x 



47T ° iX m°iXm + IXrn Xv 



= a, (-|^ TOn + ,4!), 



Therefore, all quantum corrections cancel exactly between the boson and ghost fields (a characteristic feature of 
supersymmetry), and the ensemble of dyons is basically governed by a classical field theory (35[. 

To find the ground state we examine the fields' potential energy being —Atk f ii m dT j dw m which we prefer to write 
restoring v m = fi m +i — H-m and T as 

V = 47rfJ2" m e Wm - Wm+1 - (50) 

m 

For constant fields w m , this is multiplied by the volume, therefore one has to find the stationary point for any given 
set of f m 's. It leads to the equations 

&P dV 

— — =4irf (v ie wl - W2 -v N e WN - wl ) = 0, — — = 4tt/ (v 2 e W2 - W3 - v 1 e m - W2 ) = 0, ... (51) 

OWi OW2 

whose solution is 

„ W ,- W2 {viv 2 v^...v n )tt ^- w . {v 1 v 2 v 3 ,...v n )tt 



etc. (52) 

The solution corresponds to all terms in eq. (|50[) being equal, despite a priori non-equal i^ m 's. Putting it back into 
eq. (|5D|l we obtain 

V = A t-KjN{v 1 v 2 ...v N )Ti 1 vx + v 2 + ... + v N = 1. (53) 
The maximum is achieved when all i/'s are equal: 

v x = v 2 = ... = v N = 1 P max = 4tt/. (54) 
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Equal v's correspond to the "maximal non-trivial" or "confining" holonomy, see eq. ([5]). Since there are no quantum 
corrections, the free energy of the dyons ensemble is simply proportional to the classical potential energy, F — —VV. 
Therefore, the maximum of V corresponds to the minimum of the free energy. Thus the free energy of the grand 
canonical ensemble has the minimum at the confining values of the holonomy (see the Introduction) . In the minimum 
the free energy is 



F a 



-AnfV 



16^ 2 4 V 
a i T 



N 2 A 4 V 



A 



a s N _ g 2 N 
" 8^2"' 



2tt 



(55) 



and there are no corrections to this result. In the last equation we have introduced the N- independent 't Hooft 
coupling A. 

Let us make a few comments. First, the free energy (|55|) has the correct behavior at large N. Second, V/T = V^' 
is in fact the Ad volume of the R 3 x S 1 space. Although we do not expect our theory to be valid at small temperatures 
(where the measure we use for same- kind dyons is probably incomplete), eq. (|55[) can be formally extended to the 
zero-temperature limit, as it correctly reproduces the extensive dependence on the Ad volume. Third, eq. (l55|) gives 
in fact the density of dyons. One can introduce separate fugacities f m for dyons of the m th kind into the partition 
function (|36p ; then the average number of dyons is found from the obvious relation 



<K m >= 



dlogZ 



91og/„ 



fm=f 



With separate fugacities, the result (|55|) is modified by replacing / — > (/1/2... /jv) 1/W 5 hence 



<K m >= 



f m =f 



N 



AirfV 



N A 4 
An 2 A 2 



(56) 



i.e. a finite (and equal) density of each kind of dyons in the 4- volume, meaning also the finite density of the KvBLL 
instantons. From the 3-dimensional point of view, the 'id density of dyons (and KvBLL instantons) is increasing as 
the temperature goes down: there are more and more instantons sitting on top of each other in 3d but spread over 
the compactified time direction. 

Let us add a few comments of speculative nature as they extend what is actually done here. We attempt to make 
contact with the phenomenology of the pure glue SU(N) Yang-Mills theory. In the real world there must be as 
many anti-self-dual dyons in the vacuum as there are dual ones, up to thermodynamic fluctuations ~ yV. For a 
crude estimate, we make the simplest assumption that adding anti-self-dual dyons just doubles the free energy. If the 
topological angle 9 is introduced, one has to change dyon fugacities / — ► / e %6 l N and anti-dyon fugacities / — ► / e~ lS l N , 
such that the KvBLL instanton whose fugacity is / acquires a phase e 10 and the anti-instanton acquires a phase 
e~ l9 [37l | . After minimization in w m and v m which goes as before, the free energy (|55[) becomes 



F = -An f 2 cos— V = - 
N 



16tt 2 



A 4 2cos^ 
N T 



(57) 



leading to the topological susceptibility 



}l>= I d A x 



' Tr FF(x) TrFF(O) 



1 d 2 F 



16tt 2 



16tt 2 



de 2 



0=0 



32tt 2 4 _ 1 A 4 
7VV ~ 27T 2 A 2 



(58) 



We see that the topological susceptibility is stable at large N as it is expected from the iV-counting rules. 
The free energy is related, via the trace anomaly, to the so-called gluon condensate [TJ H3| 



Fc gg l?gg± yffl 



12 



16tt 2 



from where we find 



<TrF 2 „> 
16tt 2 



N 



12 1 A 4 

n 27T 2 a 2 
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(59) 



It is the expected A^-dependence of the condensate. 

As the temperature increases, the perturbative potential energy ^ becomes increasingly important since its contri- 
bution grows as T 4 with respect to the non-perturbative one. The perturbative energy arises from the small-oscillation 
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determinant [det(— D 2 )] 1 denoted as c in eq. ([55]) . If we naively add up the dyon-induced free energy (f53|) and the 
perturbative energy (Tj) both computed at the maximally non-trivial holonomy l|54p , we obtain the full free energy 

327T 2 A 4 , T 3(2^^-l\ 



g 4 T 180 N 2 

It becomes positive and hence less favorable than the zero energy of the trivial holonomy at the temperature 

4 45 iV 4 A 4 

At this temperature, the deconfmement phase transition is expected. We see that T c is stable in N as it should be 
on general grounds. For a numerical estimate at N = 3, we take A = 1/4 compatible with the commonly assumed 
freezing of a s at the value of 0.5, and A = 200 MeV in the Pauli-Villars scheme. We then obtain from eas. (|5"51l6T)|) the 
topological susceptibility, the gluon condensate and the critical temperature (189 MeV) 4 , (255 MeV) 4 and 278 MeV, 
respectively, being in reasonable agreement with the phenomenological and lattice values. More robust quantities 
(both from the theoretical and lattice viewpoints) are those measured in units of the string tension; such comparison 
will be made in the next Section. 

From now on, we shall assume we are far enough below the critical temperature, so that the minimum of the free 
energy implies the "confining" holonomy, eq. (|54p . From eq. (|52|) we learn that at the minimum all constant parts 
of WmS are equal (up to possible difference in 2irik with integer k, which does not change the exponents of w). Let 
us note that had we imposed the overall neutrality condition of the dyon ensemble by an additional integration over 
the 9 angles (see eq. (f3Tj)) it would be equivalent to shifting w m — » w m + i9 m . Since in eq. (|46)l one integrates over 
all functions w m including their constant parts, an additional integration over #'s is unnecessary, and the neutrality 
condition is imposed automatically. 

The triviality of the free energy (|55p (which is due to the cancelation between boson and ghost quantum deter- 
minants) does not mean the triviality of the ensemble: dyons are in fact strongly correlated, as we shall see in the 
last sections. To study correlations, one has to insert source terms into the partition function (|4"6")l . With the sources 
switched on, the fields w rn are allowed to be x-dependent. Therefore, one has to retain the term —Airf^imdT/dwm 
which we rewrite using eq. (f5T)j) and eq. ([51]) as 

Action = J d 3 x ^-T{w) (61) 

where !F(w) is defined in eq. (|45|) . 

Finally, we note that the equation of motion for the fields w m , following from the (^-function (|47|) . is known as 
the periodic Toda lattice [38| which has plenty of soliton solutions. In particular, there are many one-dimensional 
domain- wall solutions interpolating between w m — w n = 2irik mn and 2nik / mn where fc, k' are integers. Why do not 
they contribute to the partition function? The answer is that any soliton is x-dependent, and an overall shift of the 
soliton is a zero mode of the operator (|48p resulting in an integration over the soliton position in space. However, it 
is also a zero mode of the identical operator for ghosts (j4"9"|) , leading to a vanishing ghost determinant. Therefore, any 
soliton gives a zero contribution to the partition function. However, solitons may and will generally contribute to the 
correlation functions. 

V. CORRELATION FUNCTION OF POLYAKOV LINES 

In the gauge where Ai(x) is chosen to be time-independent the Polyakov line is Tri(z) = Tr exp(iAj(z)/T). The 
A4 field of K KvBLL instantons away from their cores is Abelian [23| and can be gauge-chosen to be diagonal: 



A 4 (z)/T = 6„ 



1 1 



2T \ l z ^mi| l z 1. 

Comparing it with eq. (|40p we observe that A 4 can be written as 



A 4 (z)/T = diag (27T/x OT - |u m (z)) , TrL(z) = ^ exp ( 27ri/i m - -v m (z) j . (62) 
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Therefore, to compute the vacuum average of any number of Polyakov lines, one has to add a source term to the 
partition function (|46[) : 



^ exp e mi ^27ri/x mi -^ Jdxv m (x)5 mmi S(x-Zx) S j + e m2 ^2mfj, m2 -^ Jdxv m (x)5 mm2 5{x-z 2 ) 



(63) 



where Z1.2... are the points in space where Polyakov lines are placed and e mi ,m 2 ... — il depending on whether one 
takes L = exp(zA 4 /T) or = exp(-iA±/T). 

The source term is linear in v m which means that integration over v m in the partition function with a source 
produces a ^-function (|4"T|) as before whose argument is now shifted by the source: 

/ y Qjr j i \ 

Dv m — > S [ — — dw m + f- e mi -S(x-z 1 )Smm 1 - e m2 -S(x-z 2 )S m m 2 - ■ ■ ■) ■ (64) 

\ 47T Wlfm 2 2 / 

The correlation function of any number of widely separated Polyakov lines in the fundamental representation is given 
by the path integral with ^-functions: 



(Tri( Zl )TrL(z 2 ) . . .) = £ e »«H»» 1 +i* ma +~) f Dw m exp (fdx^-T(w) 

mi,m 2 ,... ^ 

n.( T 2 dT i . i . \ / T 2 9 2 JF 

M"4^ 9 Wm+f d^Z -^(x-z 1 )6 mmi --S(x-z 2 )6 mm2 -... deti--d 5 mn + f^—^ 



(65) 



It is understood that eq. (|6"5|) is divided by the same expression but without the source. The last factor comes from 
integrating over the ghost fields. 

The strategy is to find all possible solutions of the (^-functions, substitute them into the action (fBTj) and to sum 
over 777.1,2.... Note that, whatever functions w m solve the 5-functions, the Jacobian arising from those (5-functions is 
again canceled exactly by the ghost determinant. Therefore, there will be no corrections to a classical calculation. 



A. Average of a single line 

The average < TrL > is expected to be zero for the confining holonomy but let us check how it follows from the 
general eq. (f6"5)) . In this case there is only one 8- function source in eq. (|6"5|) . One has to solve the equation 

--^-d 2 w m + f^— = I <5(x-zi)<5 mmi 
47T ow m 2 

and plug the solution into the action (|6"Tj) . The solution is w m (x) w S mmi (i/2T)/\x — z±\ near the source where 
the Laplacian is the leading term and dT jdw m can be neglected. At large distances from the source w m decays, 
therefore dT /dw m can be expanded to the linear order in w m . The solution decreases exponentially with the distance. 
At intermediate distances the non-linearity is essential. However, whatever is the precise form of the solution of this 
non-linear equation, the action is finite and independent of mi, as there is a perfect cyclic symmetry in mi. Therefore, 
the action factors out from the summation over mi, and we obtain 

< TrL >= const. exp (2ni[i mi ) = 0, (66) 

mi 

as expected in the confining phase. We use here the "maximally-non-trivial" holonomy ([5]) which has been shown in 
Section IV to bring the free energy to the minimum. 

B. Heavy quark potential 

The correlation function of two Polyakov lines in the fundamental representation at spatial points zi and z 2 is 

(67) 



(TrZ(zi)Tr£t( Z2 )> = ^ e*^'-^' / Dw m exp ( I ' dx^T{w) 



mi ,71i 



II S ( ~^ d%Wm + f~^~~ ~ \ (J ( x_z i)^mi + \ S(x-z 2 )S mTll \ det (~-^d 2 6 mn + f ^ ^ 
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We are interested in the asymptotics of the correlator (|6T[) at large source separations, \z\ — Z2 1 — ► oo. We shall see 
in a moment that u> m 's solving the <5-functions fall off exponentially from the sources, w ~ exp(— M|z— z li2 |)/|z— z li2 |, 
therefore the generally non-linear equations on w m can be linearized far from the sources. The same Yukawa (or, more 
precisely, Coulomb) functions are the solutions close to the sources, as the leading term there is the Laplacian, and 
the dJ- '/dw m term can be neglected. In the intermediate range the non-linearity is essential but it has no influence 
on the asymptotics of the potential between two infinitely heavy quarks, - only on the residue of the correlator. The 
action acquires the |zi — Z2|-dependent contribution from the range of integration far away from both sources where 
w m (x) is small. Therefore, to find the asymptotics of the heavy-quark potential one can take dJ- /dw m to the linear 
order in w m and !F(w) to the quadratic order. We have for small w m 



N - 



- w m M inn IV n , 



dT 

dw m 



M 



inn u-'n 3 



(68) 



where M is the matrix made of scalar products of the simple roots of the gauge group, supplemented by a non-simple 
root to make it periodic: A4 mn = Tr C m C n> see eq. <jTJ . In our case of SU(N) 



/ 2 -1 
' -1 2 -1 



M = 0-1 

V-i '0' 

The SU (2) group is a special case where this matrix is 

= 



2 

-2 



\ 



■1 2 J 



(69) 



(70) 



The ortho-normalized eigenvectors are the pairs 



y(M) 



2 (2itk \ Trffe2l Pi . /27rfc 

— cos I n I and V} ' ' = \ j — sm ( 



V N 



N 



(71) 



corresponding to the twice degenerate eigenvalues 



M iL > = (2sin^ 



k = 1. 



JV-1 



(72) 



There is also an eig envector vj; 0) = cos(2ir-0/N-n)/VN = (1, 1, 1)/VN with a non-degenerate zero eigenvalue, and 
in case of even N there is an additional eigenvector Vn* = cos(27r- N/2/N-n)/\fN = (1, — 1, 1, —l)/*/N with a 
non degenerate eigenvalue equal 2 2 . In other words, the eigenvalues are 



2sin^ ) . /,' = 0. . 



,JV-1, 



(73) 



where the pairs of eigenvalues corresponding to k and N — k are apparently degenerate. 
In the linearized form, the (^-functions in eq. (|67p impose the equations 



d 2 w m + M 2 M m v 



2ni 

w n = — (5 mmi 6(x-zi) - 6 mni 8(x-z 2 )) 



where we have introduced the 'dual photon' mass 



(74) 



, r2 4tt/ 16tt 2 A 4 _, a ~ 



T g 4 T 2 

Equations (|T4[) are best solved in the momentum space: 

1 



2ni 

"m(Pj = — 



APM 



E„(p), 



where E„ = 5 nmi e ip ' Zl - S nni e lp ' Z2 . 



(75) 



(76) 
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This must be put into the action (|6Tj) where T{w) is to be expanded to the quadratic order. We have 



,3 _4vr/ 1 ... ... 2nf f d 3 p 



d x ~^- ■^w m {yL)M mn w n {yi) = J 3 w m (p)X m „w n (-p) 

2nf (2tt) 2 f d 3 p ( 1 \ ( 1 \ , 

= --jr -jr Jj^ e -^ { p 2 + APM ) mp M„ {^ + M 2 M ) qn M-p) 

where we have diagonalized the matrices by the orthogonal transformation built of the eigenvectors Vm corresponding 
to the eigenvalues M."'. We now pick from E m (p)E n (— p) the cross terms depending on Zi — z 2 as only they are 
relevant for the interaction. The inverse Fourier transform is 



^3p e ip-(zi-z 2 ) i e -|*i-»2 1 MVmW 



(2tt) 3 (p 2 + M 2 7W('))2 8tt MVM« 
Therefore, we continue the chain of eqs. (|77p and write 

27T 2 / 



m = ^4 V& VA^ e -l---l^^ F « (78) 

where summation over all eigenvalues labeled by I is implied. The coefficient -2 arises because cross terms in E m E n 
have negative relative sign, and there are two such terms. We obtain from eq. (|67p 

(TrL( Zl )TrLt(z 2 )) = £ exp U^^-^ + ^J^ V«\ e~^- *\**Fm V w\ . (79) 

mi.ni ^ ' 

At large separations |zi — z 2 | between the point sources the second term in eq. ([79")) is exponentially small and one 
can Taylor-expand it. The zero-order term is zero(!) as it is the product of two independent sums over m\ and n\, 
i.e. it is the product of unconnected <TrL(zi)xTrL^(z 2 )>= 0, as explained in the previous subsection. In the first 
non-zero order we get 

(Tri(zOTrLt(z 2 )) = ] ^^v / ^ e -' z --' MV ^ T7T £ exp (2^( Mmi - M „J) V/« ^ • (80) 

l— 1 mi,ni-l 

It is a sum of exponentially decaying contributions with the exponents determined by the eigenvalues see 
eq. ([75)) . The weight of the I th contribution is determined by the summation over mi,ni. For I — 1,..., [^j^l , 
eigenvalues are twice degenerate and we use the eigenvectors ([7T)) . At even TV the highest eigenvalue is non-degenerate, 
the corresponding eigenvector being Vn N ^ 2 ^ — cosirn/y/N. Summation over mi,nx in eq. (|80j) gives 

E( m\ — ni\ 2 / 2irlmi 2irlni 2itlm\ 2irln\\ 
exp 27ri — — — cos — — — cos — — — hsm — — — sm — — — (81) 
V JV / JV\ JV N N N i 

N-l 



mi , n i — l 

= NSi A , forallZ = l 

JV 



, any TV; 



5^ exp 2 



mi,ni-l 



mi — n{\ 1 

TTl C0S7TTO1 COS7Tni 

N JN 



N 

= N <5jv,2 > f° r i = "2", even TV. 

We see that only the exponent with the lowest eigenvalue y/A4^ = 2 sin contributes in eq. (|79p to the correlator of 
Polyakov lines in the fundamental representation; higher eigenvalues decouple through orthogonality. We thus obtain 

<TrL( Zl )TrLt( Z2 )) = _?l!Z_ 2 sin £ TV exp (-| Zl -z 2 | M 2 sin , (82) 
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plus exponentially small corrections from the expansion of eq. (I79|) to higher orders. This should be compared with 
the standard definition of the heavy-quark potential 



(Tri(zi)Tri t (z 2 )) = C exp 



F(zi-z 2 ) 
T 



from where we deduce the linear heavy-quark potential at large separations: 



V(z 1 -z 2 ) = |z 1 -z 2 |MT2sin- =a\z 1 -z 2 \, C = O(N ), 



(83) 



with the 'string tension' 



a = MT 2 sin — 

N 



T 



4tt/ 
T 



A 2 



A 2 N 



2 sin — = 8n sin — = — sin — 



N 



9' 



N 



A 7T 



N 



(84) 



In the last equation in the chain the A^-independent 't Hooft coupling A has been used. We see that the string tension 
turns out to be i) independent of the temperature [30( ] and it) independent of N at large N, as expected. In reality 
we expect that anti-self-dual dyons not accounted for here double M 2 and hence the dyon-induced string tension is 
actually y/2 times bigger. A more robust quantity (both from the theoretical and lattice viewpoints) is the ratio T c ^fu 
since in this ratio the poorly known parameters A and A cancel out, see eq. ()60[) : 




JV->oo 



o 



1 

TV 2 



(85) 



The values are compared to those measured in lattice simulations of the pure SU(N) gauge theories [4Q| in Table 1 
demonstrating a good agreement. The relatively large 4% deviation for the 5/7(2) group may be related to the fact 
that we have determined T c in Section IV by comparing the free energy for confining and trivial holonomy, that is 
assuming a first-order transition, whereas for N = 2 it is actually a second order one. 



N = 2 


3 


4 


6 


8 


oo 


(0.7425) 


0.6430 


0.6150 


0.5967 


0.5906 


0.5830 + «f^ + M2°« + .. . 


0.7091(36) 


0.6462(30) 


0.6344(81) 


0.6101(51) 


0.5928(107) 


0.5970(38) + ^221 (fit) 



TABLE I: Deconfinement temperature T c /yfa from eq. (|85p (upper row) and from lattice simulations (4(j (lower row) 



In Table II we add the comparison of the topological susceptibility ([58]) measured in units of the string tension, 
with the lattice data [4l|. The agreement is also remarkably good, given the approximate nature of the model. 



N = 2 


3 


4 


5 


oo 


0.5 


0.439 


0.420 


0.412 


0.399 + + M|2 + ... 


0.4831(56) 


0.434(10) 


0.387(17) 


0.387(21) 


0.376(20) + (fit) 



TABLE II: Topological susceptibility <Qr > 1//4 I \fa from eq. (I58p (upper row) and from lattice simulations [4lJ (lower row). 
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C. iV-ality and fc-strings 



All irreducible representations of the SU(N) group fall into N classes: those that appear in the direct product 
of any number of adjoint representations, and those that appear in the direct product of any number of adjoint 
representations with the irreducible representation being the rank-A antisymmetric tensor, A = 1, . . . , TV— 1. "iV-ality" 
is said to be zero in the first case and equal to A in the second. iV-ality-zero representations transform trivially under 
the center of the group Zjy; the rest acquire a phase 2wk/N, 

One expects that there is no asymptotic linear potential between static color sources in the adjoint representation 
as such sources are screened by gluons. If a representation is found in a direct product of some number of adjoint 
representations and a rank-A antisymmetric representation, the adjoint ones "cancel out" as they can be all screened 
by an appropriate number of gluons. Therefore, from the confinement viewpoint all iV-ality = k representations are 
equivalent and there are only N — 1 string tensions cr(A, N) being the coefficients in the asymptotic linear potential 

AM ' 
fel(iV-fc)! 



for sources in the antisymmetric rank-A representation. Its dimension is d(k, N) — tttztztu and the eigenvalue of the 



quadratic Casimir operator is C(k, N) — k(N — k). 

The value k = 1 corresponds to the fundamental representation whereas k = N—l corresponds to the representation 
conjugate to the fundamental [quarks and anti-quarks]. In general, the rank-(iV— k) antisymmetric representation is 
conjugate to the rank-A: one; it has the same dimension and the same string tension, u(k, N) — a(N—k, N). Therefore, 
for odd N all string tensions appear in equal pairs; for even AT, apart from pairs, there is one privileged representation 
with k = tj- which has no pair and is real. The total number of different string tensions is thus [-y] . 

The behavior of er(A, N) as function of k and N is an important issue as it discriminates between various confinement 
mechanisms. On general iV-counting grounds one can only infer that at large N and k <C N, a(k, N)/a(l, N) = 
(k/N)(l + C(l/iV 2 )) (43[. In this subsection we show that the dyon ensemble leads to the sine law for the A-strings, 

cr(k, N) = const, sin — (plus temperature dependent corrections) (86) 

satisfying the above requirement on the asymptotics. The sine behavior has been found in certain supersymmetric 
theories [44( ■ Here it follows from a direct calculation of the correlator of Polyakov lines in the rank-A antisymmetric 
representation. 

We first show that there is no asymptotic linear potential between adjoint sources. If A4 is diagonal and given by 
eq. (f62|) the eigenvalues of the Polyakov loop in the adjoint representation are exp {±i{A± m — A± n )/T), and there are 
N — l unity eigenvalues. Therefore, the average of the adjoint line is non-zero, and the correlator of two such lines 
tends asymptotically to a non-zero constant. 

Let the Polyakov line in the fundamental representation be L(z) = exp(iv44(z)/T) = diag(zi, Z2, ■ ■ ■ , Zn) where 
z rn — exp (27ri/i m — |t) m (z)J, see eq. (1621) . The Polyakov lines in the antisymmetric rank-A representation are then 

N 

L(1,N) = TrL= *m, k = l, 

m—1 

1 - 

L(2, N) = - ((TrL) 2 - TrL 2 ) = £ z m z n , A = 2, 

1 N 

L{3,N) = - ((TrL) 3 -3Tr J L 2 Tr J L + 2TrL 3 ) = ^ z mZ n z p , k = 3, 

m<n<p 

N 

L(k,N) = 2J z mi z m2 ...Zm k . (87) 

mi<m2<...<m k 

Therefore, any general L(A, N)(z) placed at the 3d point z serves as a source 



E 



7 

2lTi (fl mi + ... +MmJ- X (V mi ( z )+ ■- +Vrn k {z)) 



exp 

mK»2<...<rot 

for the v m field, which should be put into the partition function ([46]). 

To get the correlation function of two lines in fc-representation we proceed as in Subsection V. B and arrive at the 
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2tt 2 / 



generalization of eq. (|80|) : 

(Tri(fc,iV)( Zl ) Trit( fcj Ar)( Z2 )) _ 

N N 

^ exp 27rz (u mi + ... 



(j) -|a 1 -a 3 |AfVA4C) 



"MnJ 



y(0 + + y(i)l [y(J) + +V d) 

' vn. 1 1 ■ ■ • 1 r m, i. r rj i 1 ■ ■ ■ 1 * 7?. l 



«tl<«i2<...<mfc n 1 <ri2<...<nk 



In deriving eq. ([55[) it is important that the maximally non-trivial holonomy §5§ is used, leading to < TrL(k, N) >= 
0, k = 1...N — 1. Higher powers of exp (— |zi— Z2\MVMS^\ have been neglected. 

Again, the correlation function of Polyakov lines is a sum of exponentially decaying contributions with the exponents 
determined by the eigenvalues A4^ l \ see eq. (f73|) . The weight of the I th contribution is given by the sum over mi^.j 
and ni2...k- We recall the eigenvectors V™ (|7T|) and observe the following important orthogonality relation: 



E 



JV 

E 



exp 



2m 



mi+ ... +m k -ni 



-n k 



N 



v£} + ... + vV + + 



rai<m 2 <...<mi, ni <n 2 <.. . <rifc 

N Sik for all twice degenerate eigenvalues 1 = 1... [■^j^l , any jV > 2fc, 
N 5N,2k for the highest, non — degenerate eigenvalue Z = |, even jV. 



(89) 



[The orthogonality relation (|8ip is a particular case of this general one, corresponding to k = 1. The derivation of 
these relations is elementary when one presents the eigenvectors in the exponential form.] 

The above orthogonality relations imply that the correlator of the lines in rank-/c antisymmetric tensor representation 
couples only to the single exponent determined by the k th eigenvalue ^/M.^ = 2 sin N > 2k; all the rest eigenvalues 
decouple [45[ . Therefore, the correlation function ([55)1 is 



(90) 



(TrL(fc,A0( Zl ) TrLt(fc,jV)( Z2 )) = 2 sin ^ N exp (-|z x - z 2 | M 



and hence the general-fc string tension is 



trffe, N) = MTVMW = M T 2 sin — = — — sin — 
1 ' ' N A 7T N 



(91) 



as announced. Lattice simulations [4(| support this regime, whereas another lattice study [43] gives somewhat smaller 
values but within two standard deviations from the values following from cq. (|91|) . For a general discussion of the sine 
regime for fc-strings, which is favored from many viewpoints, see [431 ] . 



VI. AREA LAW FOR SPATIAL WILSON LOOPS 



The area behavior of the spatial Wilson loops is not directly related to the linear confining potential, however it is 
believed that in a confining theory the spatial Wilson loop must exhibit the area law. The reason is that i) at T — > 
Lorentz symmetry is restored, therefore the spatial loop must behave in the same way as the time-like one whose 
area law is related to the linear confining potential, ii) at high T the spatial loop eventually becomes a time-like loop 
from the 2 + 1 dimensions' point of view, which has to obey the area law to fulfill confinement in 3d. Therefore, it 
is very plausible that the spatial Wilson loop has the area behavior at any temperatures. It is expected that the 
spatial string tension is roughly constant below the deconfinement transition, and eventually grows as ~ T 2 at very 
high temperatures where the theory is basically 3-dimensional. 

In this section we demonstrate that the dyon ensemble induces the area law for spatial Wilson loops and that the 
string tension coincides with that found in the previous section from the correlators of the Polyakov lines. We think 
that it is an interesting result since a) the way we derive the string tension for spatial loops is very different from that 
for Polyakov lines, b) in a sense it demonstrates that our ensemble restores Lorentz symmetry at low temperatures, 
despite its 3-dimensional formulation. 

The condition that A4 is time-independent only partially fixes the gauge: one can still perform time-independent 
gauge transformations. This freedom can be used to make diagonal (i.e. Abelian). This necessarily implies Dirac 
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string singularities which are pure gauge artifacts as they do not carry any energy. Moreover, the Dirac strings' 
directions are also subject to the freedom of the gauge choice. In Refs. 0, [23[ the gauge choice in the explicit 
construction of the KvBLL instanton was such that the Dirac strings were connecting individual dyon constituents of 
the instanton. This choice is, however, not convenient in the ensemble of dyons as dyons have to loose their "memory" 
to what particular instanton they belong to. The natural gauge is where all Dirac strings of all dyons are directed to 
infinity along some axis, e.g. along the z axis. The dyons' field in this gauge is given explicitly in Ref. [42l ] (for the 
SU{2) group). 

In this gauge, the magnetic field of dyons beyond their cores is also Abelian and is a superposition of the Abelian 
fields of individual dyons. For large Wilson loops we are interested in, the field of a large number of dyons contribute 
as they have a slowly decreasing l/|x— x,-| asymptotics, hence the use of the field outside the cores is justified. Owing 
to self-duality, 

[A(x)] mn = [diA A ( X )] mn = -~ S mn <^ m (x), (92) 
cf. eq. (|40|) . Since Ai is Abelian beyond the cores, one can use the Stokes theorem for the spatial Wilson loop: 

W = Tr P exp i j> A i dx i = Trexp i J B l dV = exp ^-i^ J dV diV^j . (93) 

Eq. may look contradictory as we first use Bi — curlAj and then Bi — diA±. Actually there is no contradiction 
as the last equation is true up to Dirac string singularities which carry away the magnetic flux. If the Dirac string 
pierces the surface spanning the loop it gives a quantized contribution exp(27ri -integer) = 1; one can also use the 
gauge freedom to direct Dirac strings parallel to the loop surface in which case there is no contribution from the Dirac 
strings at all. 

Let us take a flat Wilson loop lying in the (xy) plane at z — 0. Then eq. |93|) is continued as 



W = ^2exp(-i— [ d 3 xd z v m S(zU =^exp(i^ / 

V Z Jx,y£Aren J _ V Z Jx, 



_ d 3 xv m d z S(z) ) . (94) 

fx,y£Avc& / m V ^ Jx,y£Arc& 



It means that the average of the Wilson loop in the dyons ensemble is given by the partition function (|46[) with the 
source 



^ exp Jd 3 x v m —jj^- Oft: V € Area) 



where 0(x 1 y £ Area) is a step function equal to unity if x, y belong to the area inside the loop and equal to zero 
otherwise. Again, the source shifts the argument of the (5-function arising from the integration over the v m variables, 
and the average Wilson loop in the fundamental representation is given by the equation 



w = E / Dw ™ ex p (J d * ^f^ 1 



(95) 



nt / T o2 , dT IT d5{z) n . A, / T n2r , d z T 

[[S [-— d 2 w m + f- + — 8 mmi — ^ 6(x, y € Area) det -— d 2 6 mn + f 



1 III 1 J r\ 1 " III III 1 , " \ w ? y ^- I I . " " jriri > J ^ 

47T ow m 2 dz J \ 47T aw m ow r , 

Therefore, one has to solve the non-linear equations on u> m 's with a source along the surface of the loop, 

- d 2 w m + M 2 (e^-^+i - e^m-i-™™) = -2m 5 mmi ^ 9(x, y G Area), M 2 = *%L, (96) 

dz T 

for all mi, plug it into the action (47r//7V)J r (w), and sum over mi. In order to evaluate the average of 
the Wilson loop in a general antisymmetric rank-fc representation, one has to take the source in eq. (|96p as 
— 2iri 5'(z) (5 mmi +. . . + S mmk ) and sum over mi < . . . < m& from 1 to TV, see eq. (|87|) . Again, the ghost determinant 
cancels exactly the Jacobian from the fluctuations of w m about the solution, therefore the classical-field calculation 
is exact. 

Contrary to the case of the Polyakov lines, one cannot, generally speaking, linearize eq. (|96|) in w m but has to solve 
the non- linear equations as they are. With no source in the r.h.s., eq. |96|) is known as t he p eriodic Toda lattice and 
it is integrable for any N. It has an hierarchy of soliton solutions constructed in Ref. [H,[4^. Below we modify those 
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solutions in such a way that they satisfy eq. ((96)) with a source in the r.h.s. We call them "pinned solitons"; their 
action determine the string tensions. We obtain below for the 'magnetic' fc-string tension 

Trk A 2 N irk . , 

a(k,N)=MT 2 sin— = — -sm— , (97) 

which coincides exactly with the 'electric' string tension (|9ip found from the correlators of the Polyakov lines. 



A. Construction of string solitons 

Let us find the pinned solitons corresponding to Wilson loops in a general antisymmetric rank-fc representation of 
the SU(N) gauge group. First of all we rewrite eq. for the difference fields w m , m +i = w m — w m+ i as only the 
differences enter the action: 

-w'{ 2 + M 2 (2e w ™ - e w ^ - e w ™) = -2niS'(z) [(S hmi + . . . + <5i, m J - (S 2 , mi +- ■ .+fe, m J] , (98) 
+ M 2 (2e w ™ - e w ™ - e w ^) = -2m6'(z) [(6 2 , mi +. . .+<5 2>m J - (<J 3 , mi +• • , 

+ M 2 (2e U) "" 1 '~ - e -«-=.«-i - e «w.i) = -2mS'{z) [(^_i, roi + . . .+<5jv-i,mJ - [S N , mi +. . .+S N , mh )} , 
-w% A + M 2 (2e w ^ - e w ^ N - e w ™) = -2ni6'{z) [(6 N , mi +. . .+5 N , mk ) - (6 hmi + . . . + S hm J} . 

At all z except one point z = where there is a source, w m .m+-\ satisfy free (zero source) equations. Solutions for very 
similar equations have been constructed in Ref. [H, [48[. Adjusting them to our case, we write the general solutions 
of eq. (f9"5| with the zero r.h.s.: 

[1 + 7 ^(m-l) £ W (z) ] [1 + 7X fc(m+l) £ W (z) ] 
W rn 771 + 1 \ Z ) = m 9 ! « = 1, ...,iV — 1, (99) 

EW(z) = exp(-MVMW z), VM^ = 2sin^, x = exp f±^) • 

The complex parameter 7 is arbitrary; there is also freedom in choosing the sign of the phase of x. Equations (|99p 
describe free "domain wall" solitons on the axis — 00 < z < +00; it is easy to check that their actions do not depend 
on 7 and are given by eq. (|9T[) . Were 7 real it could be understood as an overall shift of the domain wall by zq where 
(zq) — 7. However, we need pinned soliton solutions with a non-zero S'(z) source in the r.h.s. To find these, we 
have to take one free solution at z > and another free solution at z < 0, where the two solutions may differ only by 
the value of the so far arbitrary complex parameter 7, and by the sign of the phase of xr. 

We fix x: = exp at z > 0, and xr = exp (— ^jp) at z < 0. The condition that only the imaginary parts of 

w m ,m+i have jumps at z = requires that 7 is a pure phase factor, 7 = e la at z > 0, and 7 = er la at z < 0. Indeed, 
with these definitions, w^ m+1 at z < are complex conjugates of the same functions at z > 0; the real parts of 

(k) n 

w m m+i are continuous functions at z = whereas the imaginary parts may have jumps. We define the logarithms 
in eq. |99|) such that they have cuts along the negative axis. Thus, the general form of a "pinned" soliton solving 
eq. (|9"5j) with a non-zero source is 

[i +7X *(».->) e W (2) 1[ 1+7X m™+i) £ W (z) ] 

m f : ,, ; , , Z > U, 

[1V"E«W] 

<Ui( z )= m=l,...,N, k = l,...,N-l. 

111 ^I^iw^ ' z < °' 

(100) 

We have now to choose the phase factor 7 such that the functions w^ m+1 {z) have ±27ri jumps at z = in accordance 
with the source in the r.h.s. of eq. (|98|) . We note that at z — > ±00 the arguments of the logarithms tend to e 1 ' , hence 
all functions tend to zero at z — > ±00. As one varies |z| from 00 to 0, the arguments of the logarithms travel in the 
complex plane, ending up at the real axis at z — > 0. The trajectories for z > and for z < are mirror images of 

one another since w^\ n+1 (— z) — fio^ m+1 (z)J . If at z — > the trajectories end up at the positive semi-axis, the 

function has no jump since the logarithm is uniquely defined there. If at z — > the trajectories end up at the negative 
semi-axis, the function has a ±27ri jump owing to the cut of the logarithm along the negative semi-axis. The sign of 
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the jump depends on whether the trajectory approaches the cut from above or from below. For given N and k, the 
only handle ruling the behavior of the trajectories in the complex plane is the phase factor 7. We shall show below 
that one can find 7 such that a given function w m , m +i has a needed jump. But before presenting explicit solutions 
for k = 1,2,... let us show that the string tension for a general k representation is given by eq. (|57|) . 
To find the string tension, one needs to compute the action on the solution (|100[) : 

a f r N 

Action(fc, N) = / d 3 x ^ [exp ^^ > ) TO+1 ) - 1 9{x, y £ Area) = -u(k, N) Area 

J m=l 

where we have subtracted the constant related to the vacuum. It is understood that the solution (|100|) is valid for 
x, y inside the loop because of the ^-function in that we have omitted for brevity; outside the loop there is no 
source and w m , m+ i{z) = is compatible with the equation. At the loop boundary w m , m +i(z) interpolates between 
(|100jl and zero. Substituting the solution (|100|) we obtain for the string tension 



,<*,*, . r dz ± [exp (04-.] -%±r d > (*'- 2+ -') 7. 

J-co m _i L ™— 1 J-00 



m—1 m — 1 



where for z < one has to change E^ k \z) — > E^ k \~ z) = 1/E^(z), x — ► x* = 1/x, 7 — * 7* = 1/7. In fact the 
integrand is invariant under such change, therefore one can proceed with the above expression integrating from —00 
to +00: the integral is equal 1/(M a/M^) and does not depend either on m or 7. Therefore there are N equal terms 
in the sum and we obtain finally the string tension 



nk 

21/ sin $ 



cr(fc, A) = 4tt/ — , \ k = AfT 2 sin 



as announced. 



B. Wilson loop in the fundamental representation, k — 1 

It is easy to verify that if one takes arg(7) inside one of the N equal-length intervals covering the whole 2ir range, 
, . / 2vr\ / 2tt 4tt\ / 4tt 6tt\ / 2tt \ 

arg(7) G 7T, 7T , 7T , 7T , 7T . 7T , • • • , — 7T H , — 7T , (101) 

&y " V N J \ N N J \ N ' NJ' '\ N' )' y ' 

eq. (|100p for fc = l gives the solutions of eq. corresponding to a single source at mi = 1, 2, 3, . . . , N, respectively. 
For example, taking arg(7) inside the first interval makes the functions W12 and wn,i discontinuous at z = where 
their imaginary parts have a 2tt jump in accordance with the source term 2iri S'(z) 5 m \. All the rest functions are 
continuous. When one moves arg(7) to the second interval in (|101[) . the functions w\2 and W23 have jumps in 
accordance with the source term 2iri 5'(z) 5 m 2 while all other functions are continuous, and so on. An example of 
the solutions for N — 3 is shown in Fig. 1 where arg(7) is taken from the first interval, in this case (tt, j). When 
arg(7) is taken from the second interval (f ; — f ) or from the third interval (— f , — tt, ) the functions change cyclically 
W U -> w 23 -* w 31 -> w 12 . 



wl2 w23 w31 




FIG. 1: A bunch of profile functions W12 (left), W23 (middle) and W31 (right) inside the 5(7(3) string for five values of the 
parameter 7: arg(7) = (4,5,6,7,8) x (vr/9). The red solid curves display imaginary parts and the blue dashed curves display 
real parts of 1012, W23, W31, respectively, as functions of the distance z from the Wilson loop plane. The string tension (the 
action) is identical for all five curves. 
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The action density also varies as function of arg(7) but is periodic with a period of At N points in the middle 

of the intervals ()101|) . namely at arg(7) — tt — w \7 , the action density is real; otherwise it is generally complex. 
It is remarkable that the action itself, or the string tension, is real and does not depend on 7. It means that arg(7) 
is a new string Goldstone mode, if one allows 7 to be a function of 2c? string coordinates - in addition to the usual 
Goldstone modes associated with long-wave deformations of the string surface. 

C. Strings for higher representations, k > 2 

Wilson loop in the antisymmetric rank-fc tensor representation is a source for fc functions w mi > ■ ■ ■ , w mk where the 
numbers mi < . . . < m^ can lie anywhere on the circle (1,2,..., N). However, Toda equations have solutions not for 
all configurations of mi . . . m^. Configurations with no classical solutions presumably give much smaller contributions 
to the Wilson loop at large areas than configurations that do generate solitons as they are stationary points. 

The strategy for finding pinned solitons corresponding to Wilson loops in higher representations is simple: one 
takes the general solution (|100[) at certain value of fc and varies the phase of 7 from tt to — tt. For any k there will 
be continuous intervals of arg(7) for which the functions w^ m+1 (z) satisfy Toda equations with a S'(z) source 
in the r.h.s. corresponding to certain sets of numbers m\ < rri2 < . . . < m^. For all intervals of arg(7), the pinned 
soliton action and hence the fc-string tension is given by eq. ([97]) and is thus degenerate in arg(7). 

We did not attempt to enumerate systematically the rapidly growing variety of solitons at arbitrary TV and k. We 
find it more instructive to describe all pinned solitons of the ST/ (6) group which is sufficiently "rich" as it possesses 
non-trivial strings with k — 1, 2 and 3. 

For k = 1 the solutions have been in fact given in subsection B: six equal-length intervals of 
arg(7) € (tt, ^) , (^-, §) , (f, 0) , (0, -§) , (-§, — ^) , (— Tp, — tt) correspond to solutions with a single source 
placed at mi = 1, 2, 3, 4, 5, 6, respectively. In all cases the string tension is tr(6, 1) = 2MT sin ^. 

For k = 2 three equal-length intervals arg(7) e (tt, -|) , (f ,— f) , (— § ,— tt) correspond to the double k = 2 sources 
at mi = l,m2=4; mi = 2,m2 = 5; mi =3, 1112 — 6, respectively. In all cases the string tension is <r(6, 2) = 2MT sin 2?. 

For k = 3 two equal-length intervals arg(7) S (tt, 0) , (0, — tt) correspond to the triple k = 3 sources at 
mi = l,m2 = 3,m,3 = 5 and mi=2,m2 = 4, 7713 = 6, respectively. In all cases the string tension is c(6,3) = 2MT sin^-. 

As a matter of fact k — 3, N — 6 is a particular case of the general rank-fc representation of the SU (2k) group. 
[Another example is the fc = 1 representation of the SU (2) group, which simultaneously is a particular case of a 
fundamental representation considered in the previous subsection.] For all k and N = 2k there are pinned solitons 
generated by fc sources placed at mi = 1, m2 = 3, 7713 =5, . . . if arg(7) S (tt, 0), and placed at mi = 2, = 4, 1713 = 6, . . . 
if arg(7) £ (0, — tt). The string tension is given by eq. (|97|) where one puts N = 2k, and is degenerate in arg(7). 

To summarize this section, we have shown that to find the spatial Wilson loop averaged over the ensemble of dyons, 
one needs to solve a chain of Toda equations with a S'(z) source in the r.h.s. We have solved those equations for 
any N and Wilson loop representation fc, finding pinned solitons in the transverse direction to the loop surface. The 
solutions generalize the famous double- layer solutions for the string in the 3d Georgi-Glashow model by Polyakov [34[ • 
The resulting 'magnetic' string tension is proportional to sin ^ and coincides exactly with the 'electric' string tension 
found in Section V from the correlators of the Polyakov lines. We have observed that the Toda equations with a 
source allow a continuous set of solutions for the string profile, characterized by a phase arg(7) G (— tt, tt), all with 
the same string tension. It means that in addition to the usual Goldstone modes related to deformations of the string 
surface, there must be an extra Goldstone mode related to the string profile. Therefore, the string theory is more 
complicated than given by the standard Nambu-Goto action, which may have important implications both for theory 
and phenomenology. 
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VII. SUMMARY 

Generalizing previous work on the subject, we have written down the metric of the moduli space for an arbitrary 
number of TV kinds of dyons in the pure SU (N) gauge theory. Assuming that it is mainly the metric and not 
the fluctuation determinant about dyons that defines the ensemble of interacting dyons, we have presented the grand 
partition function of the ensemble (where the number of particles is not fixed beforehand but found from the minimum 
of the free energy at given temperature) as a path integral over N—l Abelian electric potentials v m and their duals 
w m , as well as over N—l ghost fields Xm'Xm- The resulting quantum field theory of those fields turns out to be 
exactly solvable owing to the cancelation between boson and ghost loops - a feature similar to that observed in 
supersymmetric theories. It enables one to make exact statements about the dyon ensemble: to find its free energy 
and correlation functions. 

The free energy appears to have the minimum at the "maximal non-trivial" holonomy corresponding to the confining 
zero value of the average Polyakov line. Calculating the correlation functions of Polyakov lines in various iV-ality = k 
representations (where k = 1, . . . N—l) we find the asymptotic linear confining potential with the fc-string tension 
proportional to sin ^ , the coefficient being calculated through the Yang-Mills scale parameter A and the 't Hooft 
coupling A. The actual value of A has to be determined self-consistently at the 2-loop level not considered here. 
Taking A = j compatible with phenomenology we observe a reasonable agreement of the estimated deconfinement 
temperature T c , the string tension a, the gluon condensate and the topological susceptibility with what is known from 
lattice simulations and phenomenology. A more robust ratio T c / ^/a independent of A and A is in surprisingly good 
agreement with the lattice data taken at N — 3, 4, 6 and 8, given the approximate nature of the model. 

We have also calculated the string tension from the area law for the average of spatial Wilson loops for any N 
and k in our dyon ensemble. The spatial ('magnetic') string tension coincides with the 'electric' string tension found 
from Polyakov lines for all N and k. We find this coincidence interesting as it indicates the restoration of Lorentz 
symmetry at low temperatures. Since the formalism used is 3-dimensional at finite temperatures, the restoration of 
Lorentz symmetry at T— >0 is by no means trivial. 

We do not pretend we have answered all the questions and obtained a realistic confining theory as we have ignored 
essential ingredients of the full Yang-Mills theory, enumerated in the Introduction. Our aim was to demonstrate 
that the integration measure over dyons has a drastic, probably a decisive effect on the ensemble of dyons, that 
the ensemble can be mathematically described by an exactly solvable field theory in three dimensions, and that the 
resulting semiclassical vacuum built of dyons has many features expected for the confining pure Yang-Mills theory. 
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